Segmentation of salt solution and hydrated salt crystals in porous limestone

This paper presents how simultaneous phase-and-amplitude retrieval during X-ray tomographic reconstruction allows a clear segmentation of sodium sulfate solution and hydrated sodium sulfate crystals in Savonnières limestone without the need for a dopant. This technique could be very useful to study salt crystallization processes at the pore scale, in order to improve the understanding and remediation of salt damage in porous media

Crystallization of hydrated and anhydrous salts in porous limestone resolved by synchrotron X-ray microtomography

ACLInternational audienceThe crystallization processes of two anhydrous salts (NaCl and Na 2SO4) and one hydrated (sodium sulfate) salt in the pore space of a natural building stone, Savonnières limestone, are studied. We imaged the salt solution distribution before and after crystallization and the solid crystal distribution in between repeated crystallization cycles using synchrotron X-ray microtomography. This technique proves to be very useful to study salt crystallization processes at the pore scale. The use of simultaneous phase-and-amplitude retrieval during X-ray tomographic reconstruction allows a clear segmentation of sodium sulfate solution and hydrated sodium sulfate crystals without the need for a dopant. Salt crystals can precipitate under unconfined as well as confined conditions in the multiple pore systems of Savonnières limestone, depending on their interconnection. Salt solution and salt crystals are located in mechanically weak zones of the limestone, which can be linked to damage patterns observed in this stone after repeated salt weathering cycles. The distribution and the process of pore filling by salt crystals that are revealed here advance the understanding of salt damage in porous media and may open ways to perform remediation. © 2014 Elsevier B.V. All rights reserved

Morpho+ : a software package for the three-dimensional analysis of X-ray computed tomography data

Analysis software packages are a useful tool to extract quantitative information from datasets obtained with X-ray tomography. At the Centre for X-ray Tomography of the Ghent University (UGCT) (www.ugct.ugent.be) a flexible 3D analysis software package, Morpho+, has been developed. Morpho+ provides the possibility to perform 3D measurements on datasets from various research fields. The porosity or volume fraction of the volume investigated can be determined. Several segmentation algorithms are included, to subdivide the volume into different objects such as pores or grains. For these objects various parameters can be extracted, such as the surface, diameter of the maximum inscribed sphere, diameter of the minimum circumscribed sphere, equivalent diameter, mean gray value, orientation and a measure for the sphericity. Additionally, it is possible to determine the connectivity by means of the Euler Number, size of the bottlenecks and the number of neighbouring objects. A lot of effort has been made to develop new algorithms (such as a new skeletonization algorithm and a method to reduce oversegmentation) and to improve existing algorithms, to ensure that large datasets can be processed in a small amount of time. The user interface is intuitive and flexible. Every step in the analysis process is visualized to enable easy interpretation of the result

Validation of in situ applicable measuring techniques for analysis of the water adsorption by stone

As the water adsorbing behaviour (WAB) of stone is a key factor for most degradation processes, its analysis is a decisive aspect when monitoring deterioration and past conservation treatments, or when selecting a proper conservation treatment. In this study the performance of various non-destructive methods for measuring the WAB are compared, with the focus on the effect of the variable factors of the methods caused by their specific design. The methods under study are the contact-sponge method (CSM), the Karsten tube (KT) and the Mirowski pipe (MIR). Their performance is compared with the standardized capillary rise method (CR) and the results are analysed in relation to the open porosity of different lithotypes. Furthermore the effect of practical encumbrances which could limit the application of these methods was valuated. It was found that KT and CSM have complementary fields of investigation, where CSM is capable of measuring the initial water uptake of less porous materials with a high precision, while KT was found commodious for measuring longer contact times for more porous lithotypes. MIR showed too many discommodities, leading to unreliable results. To adequately compare the results of the different methods, the size of the contact area appears to be the most influential factor, whereas the contact material and pressure on the surface do not indicate a significant influence on the results. The study of these factors is currently being extended by visualization of the water adsorption process via X-ray and neutron radiography in combination with physico-mathematical models describing the WAB

3D image analysis of a volcanic deposit

During the last decades, X-ray micro CT has become a well established technique for non-destructive testing in a wide variety of research fields. Using a series of X-ray transmission images of the sample at different projection angles, a stack of 2D cross-sections is reconstructed, resulting in a 3D volume representing the X-ray attenuation coefficients of the sample. Since the attenuation coefficient of a material depends on its density and atomic number, this volume provides valuable information about the internal structure and composition of the sample. Although much qualitative information can be derived directly from this 3D volume, researchers usually require more quantitative results to be able to provide a full characterization of the sample under investigation. This type of information needs to be retrieved using specialized image processing software. For most samples, it is imperative that this processing is performed on the 3D volume as a whole, since a sequence of 2D cross sections usually forms an inadequate approximation of the actual structure. The complete processing of a volume consists of three sequential steps. First, the volume is segmented into a set of objects. What these objects represent depends on what property of the sample needs to be analysed. The objects can be for instance concavities, dense inclusions or the matrix of the sample. When dealing with noisy data, it might be necessary to filter the data before applying the segmentation. The second step is the separation of connected objects into a set of smaller objects. This is necessary when objects appear to be connected because of the limited resolution and contrast of the scan. Separation can also be useful when the sample contains a network structure and one wants to study the individual cells of the network. The third and last step consists of the actual analysis of the various objects to derive the different parameters of interest. While some parameters require extensive calculations, others can be obtained easily. The different parameters which can be obtained are related to the size, shape and orientation of the objects. Additionally, the connectivity of a network can be analysed by comparing the set of objects before and after separation. The size of each object can be characterized by its volume, equivalent diameter and the diameter of the maximum inscribed sphere. The surface can be determined by extracting a polygonal mesh from the volume data. Calculation of Feret’s diameter reveals information about the objects elongation. Additionally, the moments of inertia can be calculated to obtain the axes of an equivalent ellipsoid. This data can be used to determine the main axis and therefore the orientation of the object within the sample. Feret’s diameter and the equivalent ellipsoid are representative for the basic shape of the object. Additionally, using a routine that fills concave regions, the convex hull of an object can be retrieved to quantify the convexity. Different ratios can be defined, which compare the surface area with the volume of the object (sphericity) or the volume of the convex hull. These ratios and the convexity characterize the objects roughness and shape. The described parameters are used to characterize volcanic deposits found in the area west of Lac Pavin (lake in Auvergne, France). The samples are taken from the most recent ‘red scoria’ layer, which is believed to be the result of the latest eruption in Western-Europe. There is however, ambiguity on the origin of the layer in terms of age and placement. The aim is to fingerprint this layer in such a way that the various eruptions in the area can be distinguished from one another. Measurements of the vesicle density, volume and connectivity of the investigated deposits provide information about the intensity of the eruption. Additionally, vesicle geometry can be related to the magmatic permeability, which is essential to the dynamics of the eruption

Latest developments in 3D analysis of geomaterials by Morpho+

At the Centre for X-ray Tomography of the Ghent University (Belgium) (www.ugct.ugent.be) besides hardware development for high-resolution X-ray CT scanners, a lot of progress is being made in the field of 3D analysis of the scanned samples. Morpho+ is a flexible 3D analysis software which provides the necessary petrophysical parameters of the scanned samples in 3D. Although Morpho+ was originally designed to provide any kind of 3D parameter, it contains some specific features especially designed for the analysis of geomaterial properties like porosity, partial porosity, pore-size distribution, grain size, grain orientation and surface determination. Additionally, the results of the 3D analysis can be visualized which enables to understand and interpret the analysis results in a straightforward way. The complementarities between high-quality X-ray CT images and flexible 3D software are opening up new gateways in the study of geomaterials